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Abstract 
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(N 

<— I , In this paper we study the well-posedness and long-time dynamics of a diffuse-interface 

model for the mixture of two viscous incompressible Newtonian fluids with thermo-induced 
Marangoni effects. The governing system consists of the Navier-Stokes equations coupled 
with phase-field and energy transport equations. We first derive an energy inequality that 
illustrates the dissipative nature of the system under the assumption that the initial temper- 
ature variation is properly small. Then we establish the existence of weak/strong solutions 
and discuss the long-time behavior as well as the stability of the system. 

Keywords. Phase-field model, Navier-Stokes equations, Marangoni effects, well-posedness, 
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1 Introduction 

> 

The study of interface dynamics is of great importance in the hydrodynamic theory of complex 
fluids. In the classical approaches (e.g., the sharp-interface model), the interface is usually 
■ considered to be a n — 1 dimensional free surface of zero width that evolves in time with the 

fluid. The resulting hydrodynamic system describing the mixture of two immiscible Newtonian 
fluids with a free interface usually consists of Navier-Stokes equations in each fluid domain 
with kinematic and force balance boundary conditions on the interface. On the other hand, 
the so-called diffuse- interface model (or phase-field model), recognizes micro-scale mixing of 
the macroscopically immiscible fluids and the interface represents a thin region with a steep 
transition property between two fluids (cf. Anderson et al [2j). Within this region, the fluid is 
mixed and has to store certain "mixing energy". The diffuse-interface model can be viewed as 
a physically motivated level-set method that describes the interface by a proper mixing energy. 
Compared with the sharp-interface model, the diffuse-interface model can describe topological 
transitions of interfaces (like pinchoff and reconnection) in a natural way (cf. Lowengrub et al 
[21] ) and it has many advantages in numerical simulations of the interfacial motion (cf. [8|ll8}l33] 
and references therein). 

The Marangoni effect was initially observed by Thomson [30J during the study of the inter- 
esting phenomenon "tears of wine" . Afterwards this phenomenon was defined in more detail in 
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Marangoni [22] in terms of surface tension gradients and named after the author. The Marangoni 
effect is a phenomenon whereby mass transfer occurs due to differences in surface tension. Such 
differences can either be attributed to non-uniform distributions of surfactants (cf. Mendes- 
Tatsis and Agble |23j) or to the existence of temperature gradient in the neighborhood of the 
interface (cf. Sterling and Scriven [26]). The latter is called the thermo-capillary convection or 
the Marangoni-Benard convection, which becomes more and more important in the application 
of complex fluids, liquid-gas systems and ocean-geophysical dynamics (cf. e.g., [31I11II51IIS1IM! ) • 
The conventional Marangoni-Benard convection for the mixture of two Newtonian flows can 
be described by a sharp-interface model involving the Boussinesq approximation (cf. e.g., Liu 
et al [19]) 

p(ut + (u ■ V)m) + Vp - vAu = -pegj, (1-1) 
V-u = 0, (1.2) 
9 t + u ■ V6> = kA6, (1.3) 

where u, p and 6 stand for the fluid velocity, the pressure, and the relative temperature (with 
respect to the reference background temperature 6b, which is assumed to be a constant for 
the sake of simplicity), respectively, p is the density of fluid mixture, v is the viscosity, g is 
the gravitational acceleration, j is the upward direction and k > is the thermal diffusion 
constant. We assume that the temperature-dependent density pg is described by the Boussinesq 
approximation 

Pe = p(l-a6), (1.4) 

where a is the coefficient of thermal expansion. The background density p is assumed to be a 
constant and the difference between the actual density and p only contributes to the buoyancy 
force. Interface conditions are given by 

l t + u-Vl = Q, (1.5) 
[T] • n = -aHn + (r • V<t)t, (1.6) 



where I stands for the interface length of the mixture. The kinematic condition (jl.5p indicates 
that the surface (I = 0) evolves with the fluid. Equation (|1.6p is the balance of forces on the 
interface, where H is the curvature of the interface, [T] is the jump of the stress across the 
interface, a is the surface tension, r is the tangential direction on the interface and n is the 
normal direction. 

In this paper, we shall investigate a diffuse-interface model, which was used to describe the 
thermo-induced Marangoni effects in the mixture of two incompressible Newtonian fluids. A 
phase-field variable (j) is introduced as the volume fraction to demarcate the two species and to 
indicate the location of the interface. The region {x : <p{x,t) = 1} is occupied by fluid 1 while 
{x : 4>(x,t) = —1} is occupied by fluid 2. The interface is represented by {x : (f>(x,t) = 0}, with 
a (fixed) transition layer of thickness e. In the diffuse-interface approach, one usually introduces 
an elastic (mixing) energy of Ginzburg-Landau type 



El 



l -\V4>\ 2 + F{4>) 



dx, (1.7) 



which represents the competition between the hydrophobic and hydrophilic effects of the two 
different species. The physically relevant energy density function F that represents the two 
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phases of the mixture usually has a double-well structure. A typical example of F is the so- 
called logarithmic potential (cf. Cahn and Hillard [6]) 



F{(j>) =ll {l-4> 2 ) + 72 [(l + ^)ln(l +4>) + (1- 0)ln(l-0)], 71,72 > 0. 

In practice, this (singular) potential is often replaced by a smooth double-well polynomial ap- 
proximation 

In this paper, we start with the simple case that the two components of the binary fluid have 
matched densities and the same constant viscosity v as well as the same constant heat conduc- 
tivity k. As a consequence, we consider the following system (cf. e.g., [18l[T9| l2"T] ) : 



p(ut + (u ■ V)u) + Vp — vAu 
= -V • 



\{6)V4> ® - A(0) ( \\V(f>\ 2 + F(<f>))I 



Pegj, (li 



.2 

V-u = 0, (1.9) 

V0 = 7(A0-F'(0)), (1.10) 

+ u • V0 = kA9, (1.11) 



for (x,t) G f2 x (0, +oo). Here, we assume that f2 is a bounded domain in R n (n = 2,3) with 
smooth boundary T. The usual Kronecker product is denoted by ®, i.e., (a (8) b)ij = aibj for 
a, b G M n . The system (|1.8p — ([l.lip contains the Navier-Stokes equations, an Allen-Cahn type 
equation for the phase-field function and an energy transport equation for the temperature. 
The parameter 7 represents the microscopic elastic relaxation time due to the presence of the 
microstructure of the mixture. As 7 — > 0, the internal dissipative mechanism will disappear 
and the limiting equation is equivalent to the mass transport equation for incompressible fluids 
(cf. [ISlEl]). The temperature-dependent surface tension coefficient X(9) is supposed to be 

A(0) = X (a-b9), 

where Ao > 0, a > 0, b ^ are constants. Usually Ao is assumed to be proportional to the 
interface length e (cf. Sun et al |27J). Numerical experiments have been made in the recent 
paper Sun et al [27] to illustrate the role played by thermal energy in the interfacial dynamics of 
two-phase flows due to the thermo-induced surface tension heterogeneity on the interface. Their 
results suggest that the system (|1.8p - (|l.lip (and its generalizations) turns out to be a suitable 
mathematical representation that reflects the thermo-induced Marangoni effects in the mixture 
of fluids. 

We suppose that the system (|1.8p - (|l.lip is subject to the initial conditions 

u\ t =o = u (x) with V • n = 0, <f>\ t =o = <j>o{x), 0\ t =o = 9 (x), i£fi. (1.12) 
Moreover, we assume no-slip boundary condition on the velocity u 

u(x,t) = 0, (x,t) G r X (0,+oo), (1.13) 
nonhomogeneous Dirichlet boundary condition on the phase function <f> 

4>(x, t) = -1, (x, t) G T x (0, +00), (1.14) 
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and homogeneous Dirichlet boundary condition on the temperature 6 



0(x,t) = 0, (x,t) e T X (0,+oo). 



(1.15) 



The goal of this paper is to provide a detailed mathematical theory of existence, uniqueness, 
regularity and long-time behavior of solutions to the non-isothermal Navier-Stokes-Allen-Cahn 
system fjl.8j) — (|X.15 j) . First, we prove the existence of global weak solutions in two and three 
spatial dimensions (cf. Theorem 12. ip . Next, we obtain the existence and uniqueness of a 
global strong solution in 2D (cf. Theorem 13. ip . a local strong solution in 3D (cf. Theorem 
13. 2|) and a global one provided that the viscosity v is properly large (cf. Theorem 13. 3|) . The 
long-time dynamics of the system seems to be more complicated than the uncoupled Navier- 
Stokes equations. We prove that as t — > +oo, the phase function converges to a solution of the 
stationary Allen-Cahn equation and the velocity u as well as the temperature 9 converges to zero 
(cf. Theorem 14. ip . Stability for minimizers of the elastic energy is also discussed (cf. Theorem 
I4.2p . We just remark that our results can be easily extended to the cases with more general 
Dirichlet boundary conditions for the phase function. For instance, ()1.14[) can be replaced by 

3 

4>{x,t)\r = h{x) on T x (0,+oo), with h(x) = (j)o(x)\r and h(x) € H2(T), \h(x)\ < 1 (we refer to 
Lin and Liu [T7] for a similar situation for a simplified nematic liquid crystal system). 

It is easy to verify that for the isothermal case of system (|1.8j) - (|1.15|) without the Boussinesq 
approximation term, there is a dissipative energy law 



This basic energy law reveals the underlying physics for the isothermal Navier-Stokes-Allen- 
Cahn system and it plays an important role in the study of well-posedness as well as long-time 
behavior of the system. We refer to the recent work [9|lll] for detailed mathematical analysis on 
an isothermal NSAC system for the incompressible two-phase flows, where long-time behavior 
of global solutions was analyzed within the theory of infinite-dimensional dissipative dynamical 
systems (e.g., the existence of global attractors, exponential attractors, trajectory attractors 
and convergence to single equiblibria) . 

However, in our present case, the surface tension parameter A in (|1.8|) depends on the tem- 
perature such that it is no longer a constant. Moreover, the Boussinesq approximation is also 
applied. These bring us challenges in mathematical analysis of the system. We are not able to 
derive the same dissipative energy equality as for the isothermal case. In particular, the special 
cancelation between the induced stress term in the momentum equation and the transport term 
in the phase-field equation is no longer valid (see Remark 12.21 below) . This relation is crucial to 
derive the dissipative energy law like (|1.16p (cf. |9l[T71[33]). Nevertheless, taking advantage of 
proper maximum principles for the phase function and temperature, we show that if the initial 
temperature variation is not large (bounded in terms of coefficients of the system), we can derive 
an energy inequality for the system (jl.8|) - (jl.l5|) . which reflects the dissipative nature of the flow. 

We remark that in the phase- field equation (|1.1U|) . the dynamics of the phase function <j) is 
assumed to be driven by a gradient flow of Allen-Cahn type. In order to keep the conservation of 
overall volume fraction, people usually assume that the internal dissipation is described through 
a Cahn-Hillard equation (with convection) for (f), which can be viewed as a gradient flow of 
the elastic energy in the Sobolev space H^ 1 (cf. |18p21j). The resulting system is nevertheless 




(1.16) 
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much more involved in mathematical analysis, because it contains a fourth-order differential 
operator and thus the maximum principle for the phase function <f> no longer holds (see, for 
instance, [U[5j[T0] for the isothermal Navier-Stokes-Cahn-Hilliard system). It seems that our 
results cannot be extended to this case in a straightforward way. For instance, in the derivation 
of the dissipative energy inequality (|2.3p . we rely on the L°°-estimate of (f), which will no longer 
be available due to the lack of maximum principle. This problem will be studied in our future 
work. 

The rest of the paper is organized as follows. In Section 2, we derive an energy inequality 
that guarantees the dissipative nature of the system and establish the existence of global weak 
solutions. In Section 3, we discuss existence and uniqueness of strong solutions in both 2D and 
3D. In Section 4, we study the long-time dynamics and stability property of the system. 

2 Global weak solutions 

If X is a real Hilbert space with inner product (-,-)x, then we denote the induced norm by 
|| • \\x- X' indicates the dual space of X and (-,-)x',x will denote the corresponding dual 
product. We indicate by X the vectorial space X n endowed with the product structure. For 
simplicity, the scalar product in L 2 (ft) (also L 2 (f2)) will be denoted by (•, •), and the associated 
norm by || • ||. For two n x n matrices Mi,M2, we denote Mi : M2 = trace(MiMrT). Let 
V = Cg°(n,R n ) n {v : V • v = 0}. We denote H (respectively V) the closure of V in L 2 
(respectively H 1 ). 



H and V are Hilbert spaces with norms || • || and || • ||jji> respectively. We recall the Stokes 
operator S : H 2 (0) nV^H such that Su = -Au + Vvr G H, for all u G H 2 (ft) n V. S' 1 is a 
compact linear operator on H and \\S • || is a norm on D(S) that is equivalent to the H 2 -norm. 
Then there exists a positive constant C = C(n, O), for which (cf. Temam |29j) 




In the following text, we denote by C, Ci the generic constants depending on a, b, g, k, Aq, 
a, 7, £1, e and the initial data. Special dependence will be pointed out explicitly in the text if 
necessary. 

Without loss of generality, we assume p = 1 in the remaining part of this paper. Now we 
introduce the weak formulation of the initial boundary value problem (|1.8p - (|1.15|) : 

Definition 2.1. For any T € (0, +00), the triple (u,4>,9) satisfying 



is called a weak solution of the problem (|1 .8j) - (jl.lip if the initial and boundary conditions (|1.12p - 
(|1.15p are satisfied and for a.e. t G (0, T), 



H = {u G L 2 : V • u = 0, u ■ n = on T}, V = {u G : V • u = 0}. 




"Hh? + IMIhi\r < C\\Su\\. 



(2.1) 



u G L°°(0, T; H) n L 2 {0, T; V), 

G L°°(0, T; H 1 n L°°) n L 2 (0, T; H 2 ), \<j>\ < 1, a.e. infix [0, T] 
9 G L°°(0,T;H^ n L°°) n L 2 (0,T; H 2 ), 



(ut, u)v',V + / ( u ' V) n ' y dx + v j Vu : Vvdx 





5 



= / [A(0)V0<8> V0] : Vvdx + ag / Oyvdx, VueV, 

</>t + u-V<f> = j(A(p - F'(4>)), a.e. in 0, 

6 t + u ■ V6> = kA9, a.e. in 0, 

u|t=o = uo(z), 0U=o = 0o(^), 6> U=o = Oo(x), in 

Remark 2.1. In order to derive the variational formulation for u, we use the following facts 
due to the incompressibility condition: for any v G V, it holds 



Vp ■ vdx = / j • vdx = 0, 
v. Jn 

n| V-(A(#)(i|V0| 2 + F(0))l 



vdx = I V 



\(e)(-m 2 + F( 



■ v dx = 0. 



Next, we state the result on the existence of global-in-time weak solutions: 

Theorem 2.1 (Existence of weak solutions). Suppose n = 2,3. For any initial data (uq,4>q,0o) 
G H x (H 1 ^) n L°°(n)) x (l#(fi) n L°°(0)) satisfying 



< 1, 



7Q||L° 



< 



4C?|6|V 2A: 



07^ 



(2.2) 



where C\ is a constant depending only on n,Q, the problem (|1.8p - (|1.15p admits at least one 
global weak solution such that 

u G L°°(0,+oo;H) nLf oc (0,+oo;V), 

<p G L°°(0, +00; H 1 n L°°) n L 2 OC (0, +00; # 2 ), |0| < 1, a.e. in O x [0, +00), 
G L°°(0, +00; Fq 1 n L°°) n L 2 OC (0, +00; F 2 ). 



2.1 Dissipative energy inequality 



An important feature of problem (|1.8|) - (|1.15|) is that <f> and 9 satisfy the following weak maximum 
principles, which will be useful in the derivation of the dissipative energy inequality for the 
system. 

Lemma 2.1. Suppose u G L°°(0, T; H) n L 2 (0, T; V). 7/0 G L°°(0, T; H 1 n L°°) n L 2 (0, T; F 2 ) 
is i/ie weaA; solution of the initial boundary value problem 

<tn + u ■ V0 = 7(A0 - F'(0)), a.e. m O, 

0(x,i)| r = -l, (x,t) GTx (0,T), 

0|t=o = <^o C^) ^ ^ n wi/t |0o| < 1 a - e - i- n ^1 

then \(j)(x,t)\ < 1, a.e. in Q for each t G (0, T). 

Proof. The proof is similar to that for the liquid crystal system (cf. e.g., [TlfTT]). so we omit the 
details here. □ 

Lemma 2.2. Suppose u G L°°(0,T; H) n L 2 (0, T; V). If 9 G L°°(0, T; #<} n L°°) n L 2 (0, T; F 2 ) 
zs i/ie u>ea£; solution of the initial boundary value problem 



9 t + u-V9 = kA9, a.e. in Q, 
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0|r = O, (x,t) e r x (0,T), 

e\ t=0 = e (x)eH^nL 00 , 

then \\9(t)\\ L °° < ||6» ||l°° for every t € (0,T). 

Proof. Multiplying the equation by \6\ q ~ 1 6 (q > 1), integrating over f2, we get 

1 d [ \6\ l +Hx+— f u-v\e\i +1 dx+ ^ q ~]l [ k v(\e\^e) 

Jn 1 + 1 Jn (1 + IF Jn 



dx = 0, 



1 + q dt Jn 
which implies that 

l|0(t)Hw < INb < ini^ll^olU-, Vg>l, ie(0,T). 

Taking the limit q — > +oo, we arrive our conclusion. □ 

In what follows, we derive a dissipative energy inequality, which turns out to be crucial in 
the study of well-posedness as well as long-time dynamics of the problem (|1.8|) - (|1.15|) . 

Proposition 2.1 (Dissipative energy inequality). For n = 2,3, we assume that the initial 
phase function <f>§ and the initial temperature 9q satisfy the assumption (|2.2p . Then there exist 
constants (,uj > that depend only on CI and coefficients of the system such that if (u,4>,9) is 
a smooth solution to the problem ()1.8|) (|1.15[) . then the following energy inequality holds: 

^<~||V«|| 2 -aAo7l|A0-^)|| 2 -A:C||Ae|| 2 <O, Vt > 0, (2.3) 

where 

£{t) = |Kt)|| 2 + aA ||V<Xt)|| 2 + 2aAo / F(<j>(t))dx + CI|V#(t)|| 2 + oj\\9(t)f > 0. (2.4) 

Jn 

Proof. Multiplying (US) with u, (OUj) with -aA o (A0 - F' (</>)), ([TIT]) with -(A9 (( > is a 
constant to be determined later), respectively, adding them up and integrating over Cl, we have 



\jt ( l|M " 2 + aA ° l|V< ^ 2 + 2aA ° J F ^ dx + Cl1 W| 
+i/||Vn|| 2 + a 7 A ||A</> - F'(^)|| 2 + KI|A0|| 2 

(A(0)V<£<8> V0) : V-udx + a / gj ■ u dx + C / (u-V)9A9dx 

Jn Jn 



n 

+a\ / {u-V)^{A^-F'{4>))dx 
Jn 

:= J1 + J2 + J3 + Ja- 

In the following we just treat the case n = 3, while the case n = 2 is similar. Recall the 
Gagliardo-Nirenberg inequality 

||V0|| L 4 < Ci(||A0||a 110111^ + ||0|| ioo ), V0 e tf 2 , 

where Ci = Ci(n, Q) depends only on n and Cl. Combining it with the Poincare inequality and 
the Young inequality, we deduce that 

J\ + J 4 
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-a\ J u • V (^y~ + Ftyfj dx - b\ J 9Vu : (V0 ® V0) dx 



< |6|Ao||e|U-||V«||||V^||i4 

< ^||Vu|| 2 + 1^011011^^^ 

< ^iivnii 2 + wcf l bl2xl \\o\\u\\^ - F'^t + wF'mmwi 

+ 8Ctl f x2 ° \\e\\lo . 

Then by Lemmas 12.11 12.21 and the assumption (|2.2p , we obtain that 

V 



Mgjg ||„|| i „|| A ^,_j.. W ||2|| ||2 M < ^|| A0 _ F . W |p 



l^HOIII.IIFWinWIi. < 18C ?y"l' ||M||i M i 

fc CiiA/»n2 54 ■ 16 4 C i L6 C||6| 8 A§|0| 



and 



£ f l|A <"' + « ( = 

£ ~8~ I|AS|1 + i¥? 1,911 ' 

where C2 depends only on 0. As a result, 

J a + J 4 < ^||V.|| 2 + 2*p||A0 - F'(0)|| 2 + ^||A0|| 2 + ^||0|| 2 , 

with 

54 ■ 16 4 C 1 16 Q||b| 8 Ag|^| 8 + 54 ■ 8 4 e ie Cl 6 C%\b\ 8 \ 8 

Next, by the Poincare inequality 

J2 < |a|M||0|||| U || < Cp|a|M||0||||Vu|| < ^||V U || 2 + \\0\\ 2 , 

where Cp depends only on 0. For J3, we have 

j 3 = -c j u ■ v O^p-) dx + c J u • [v ■ (v# ® w)] 

= -C / Vn : (V6< <g> V6<) dx 

< Cl|Vn||||V0|| 2 4 <C 2 C||Vn||||A0||||0|| L ^ 

< ^l|Afl|| a +3£||flb||i-HV«|| a 



Taking 

fc|b| 2 A 

we infer from the above estimates that 

(jM| 2 + aA o ||V0|| 2 + 2aA o J F(<j>)dx + C|| V0| 

+^||Vu|| 2 + 7 aA ||A^ - F'(0)|| 2 + KI|A0|| 2 
< ^4||^|| 2 , (2.5) 

where 

2C 3 2C 2 |a[ 2 | g | 2 
W — ~ ^~ H • 



C 2 C4 

Multiplying (OB by 2w(9, cj — — > 0, integrating over Q,, and using the Poincare inequality, 
we obtain 

w i|0||2 = _2^||V0|| 2 < _ g^p! = -C 4 ||0|| 2 . (2.6) 
Adding (|2.5|) with (|2.6|) . we arrive at our conclusion. □ 

Remark 2.2. For the isothermal case of the system (|1.8|) - (|1.15|) without the Boussinesq approx- 
imation term, there is a special cancelation between the induced stress term in the Navier-Stokes 
equations and the convection term in the Allen-Cahn equation, which yields the dissipative en- 
ergy equality (|1.16|) . However, for the current non-isothermal system f)1.8f ) — fjl. 15[) . there exists 
an extra high- order term 

-b\ / 6Vu : (V0 <g> V0) dx 
Jn 

containing velocity, phase function and temperature that cannot be eliminated (the Boussinesq 
approximation is a lower-order term and is easier to handle). In order to overcome this difficulty, 
we introduce the smallness assumption (|2.2p and try to seek certain energy dissipative inequality 
instead. 



Remark 2.3. It is worth mentioning that the conditions in (|2.2p does not involve the interfacial 
parameter e. 



2.2 Proof of Theorem I2TTI 

The proof is based on a semi-Galerkin method (see, for instance, Lin and Liu |17j for a simplified 
nematic liquid crystal system). Let {wi} (i € N) be an orthonormal basis of V formed by the 
eigenvectors of the Stokes problem 

— Awi + VPj = XiWi, in Q, Wi\r = 0, 

with \\wi\\ = 1 and < Ai < A2 < • • • < A n < • • • with A n — > 00 as n — > +00. 

For every m G N, let V m = span{w\, W2, • • ■ ,w m }. We denote by P m : H — > V m the 
orthogonal projection. Given (u , tp , 9 ) € H x (F x (fi) n L°°(fi)) x (F^(O) n L°°(fi)) satisfying 
(|2.2|) . we consider the following approximate problem: 

(d t Um, w m ) V ',v + / (v-m ■ V)u m • v m dx + v Vu m : Vv m dx 
Jn Jn 
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= / [A(0 m )V0 m (g) V0 TO ] : Vv m dx + ag Q m j ■ v m dx, \t v m G V m , (2.7) 

d t (j>m + u m -V(j) m = j(A(j) m -F\(J) m )), a.e. in O, (2.8) 

9 t m + u m • V6> m = /cA0 m , a.e. in Q, (2.9) 

<^ m (x,t) = -l, m (x,t) =0 onT, (2.10) 

u m \t=o = P m u (x), (j) m \ t =o = <f>o(x), m \ t =o = 6 {x). (2.11) 

Indeed, we observe that all the a priori bounds derived (formally) from the energy inequality 
(|2.3p still hold for the approximate problem. If we fix u m G C([0,T];V m ), then we can find 
4>m = 4>m\^m\ and 6 m = m [u m ] solving (|2.8|) and (|2.9p (with u m = u m ), respectively. Inserting 
(j) m and 8 m into the equation (|2.7p . we can find a solution u m = T[u m ] that defines a mapping 
Urn *-> T[u m \. On account of the a priori bounds, we can easily show that T admits a fixed point 
by means of the classical Schauder's argument on (0, To), with < Tq <T. Finally, applying 
again the a priori estimates, we are allowed to conclude that the approximate solutions can be 
extended to the whole time interval [0, +oo) (cf. also Ezquerra et al Appendix]). Since the a 
priori estimates of the approximate solution are uniform in parameter m, then using a similar 
argument as in |17t Section 2], we can pass to the limit m — > +oo and complete the proof of 
Theorem I2JJ The details are omitted here. 



Corollary 2.1. Suppose n = 2,3. Under the assumptions of Theorem \2.1l the weak solution 
(u, 4>i 0) to the problem (jl.8p — f|l . 15j) satisfies 



iw)iu~<i, \\m\\L~<^^, w> , 



and the energy inequality 

E(t) + J [-\\Vu\\ 2 + a\ j\\A(f>-F'(<l))\\ 2 + k£\\AO\\ 2 )dt<£(0), Vt > 0, 
which yields the following uniform estimates: 

Ht)f + U(t)\\m + P(t)\\m < M > Vt>0, 

(\\Vu(t)\\ 2 + \\AJ>(t) - F'm)\\ 2 + \\M(t)\\ 2 )dt < M, 

where M > is a constant depending on \\uq\\, H^oIIh 1 ; II^oIIh 1 ) an d coefficients of the system. 



o 



3 Strong solutions 

In this section, we prove the existence and uniqueness of strong solutions to the problem (j 1 . 8 [) — 

Definition 3.1. For any T G (0, +oo], u G V, 4> G H 2 (Q),9 G H 2 (Q) n H%(£1), we say that 
the triple (u,<p,6) is a strong solution to the problem (jl.8p - (|1.15p . if(u,<p,6) is a weak solution 
and 

u G L°°(0, T; V) n L 2 (0, T; H 2 ), 0, 6 G L°°(0, T; F 2 ) n L 2 (0, T; H 3 ). 
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Based on the semi-Galerkin scheme in the previous section, in order to prove the existence of 
strong solutions, it suffices to derive proper uniform higher-order estimates for the approximate 
solutions and then pass to the limit m — > +00. We observe that the approximate solutions satisfy 
the same basic energy inequality and higher-order differential inequalities as smooth solutions 
of the problem (|1.8|) - (|1.15|) . Thus, for the sake of simplicity, all the calculations below will be 
carried out formally for smooth solutions. 

The main results of this section are as follows. 

Theorem 3.1 (Global strong solution in 2D). Suppose n = 2. For any uq G V, <f>Q G H 2 (Q),9o G 
(-ffg (Q)PiH 2 (Q)^ satisfying the assumption (|2.2p . the problem (|1.8p - (|1.15p admits a unique global 
strong solution such that 

u G L°°(0, +00; V) n L 2 OC (0, +00; H 2 ), 

<f> G L°°(0, +00; H 2 ) n L 2 OC (0, +00; H 3 ), 

9 G L°°(0, +00; H 2 n Hi) n L 2 OC (0, +00; H 3 ). 

Theorem 3.2 (Local strong solution in 3D). Suppose n = 3. For any uq G V, 0q £ H 2 (CI),6q G 
(Hq (Q) nH 2 (Q)) satisfying the assumption (|2.2j) . the problem f|l .8j) — fjl . 15 j) admits a unique local 
strong solution. 

Theorem 3.3 (Global strong solution in 3-D under large viscosity). Suppose n = 3. For any 
uq G V, (j)Q G H 2 (Q),6o G (Hq(Q) nH 2 (Q)^ satisfying the assumption (|2.2p . if in addition, the 
lower bound of the viscosity, i.e., v is sufficiently large (cf. (|3.12p ), then the problem (|1.8p - (jl.l5p 
admits a unique global strong solution. 

3.1 Two dimensional case 

First, we are going to derive a specific type of higher-order energy inequality in the spirit of [17] . 

Lemma 3.1. Suppose n = 2. Let (u,(p,9) be a smooth solution to the problem (jl.8p - (|1.15p . We 

introduce the quantity 

Ai(t) = \\Vu(t)\\ 2 + a\ \\A^(t) - F>^(t))\\ 2 + ??i||A^)|| 2 . (3.1) 
Then the following differential inequality holds: 

^i(t) + i/||5n|| 2 + aA O 7l|V(A^-F , (0))|| 2 +7 ?1 A ; ||VA0|| 2 < C{A\{t) + A 1 {t)), Vt > 0. (3.2) 

r/i and C are two positive constants which may depend on \\uq\\, H^oHi/ 1 * 1 1 1 1 if 1 > an d coeffi- 
cients of the system. 

Proof. We observe that — (Am, Uf) = (Su,ut), since Ut G H. Besides, due to (|1.14p . we infer that 
4>t\r = 0, then by (|1.13p and the equation ()1.10p . we see that A<j) — F'(<j))\„ = 4>t + u • V0|r = 0. 
Using equations (|1.8p and (|1.10p . we compute that 

^(l|Vn|| 2 + aA o ||A0 - F'(<p)\\ 2 ) + u\\Suf + aA o7 ||V(A^ - F'(0))|| 2 
= -(u ■ Vu, Su) + 6A (V • (0V(f> <S> V0), Su) + ag(9j, Su) 

-aA o (V<A(A0 - F'(<t>)), Vtt) - a\ 0l (F"(ct>)(Act> - F' {4>)) , A<j> - F' ((/>)) 
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-2aA / {A(f>- F'i^VjUiS/jV^dx 
Jn 

6 

K m- (3-3) 



m=l 



In what follows, we proceed to estimate the right-hand side of (|3.3[) by using the uniform 
estimates obtained in Corollary 12.11 and properties of the Stokes operator. 



Kt < ^||Su|| 2 + C|M|Lo||Vu|| 2 

< ^\\Suf + C\\n\\\\Au\\\\Vnf 

< ^-\\Su\\ 2 + C\\Vu\\\ 



For the second term K 2 , we have 

K 2 = b\ [ VjOVicfNj^Su^dx + ^(#V|V0| 2 , Su) + 6A (#A</>V^, Su) 
Jn 2 

= 6A / VjOVicjNj^Su^dx - ^(|V0| 2 V0, Su) + 6A O (#A0V^, Su) 
Jn 2 

:= K 2a + K 2b + K 2c . 

K 2a + K 2b 

< C||V0|| L6 ||V0|| 2 6 || < Sn|| 

< ^\\Su\\ 2 + C(||VA0||1||V0||I + ||W|| 2 )(||VA0||§||V0||§ + ||V</>|| 4 ) 

< ^H^ll 2 -K C(||X7^||f HV^nf -k ||V6»|| 2 )(||V(^<^ - ^'(^))||f -h ||^"(^)||t_ -K C) 

< ^II^H 2 + ^||V(A^ - F'(0))|| 2 + C 5 ||VA0|| 2 + C(||W|| 2 + ||W|| 6 ), 

where C5 is a constant depending on ||«o||> H0o||.H"i> ll^ollifij ^ an d coefficients of the system. 

K 2c < |6|Ao||e|| r oo||A0||||V0||L«||5«|| 

< C||0||,H|A0||||V0||| 2 ||V#||S«|| 

< ^ll^ll 2 + <7||0||£»(||A* - F'(<p)f + C)(||V(A0 - F'm\ + II A* - F'(0)|| + C) 

< ^ll^ll 2 + C||A0 - F'(0)|| 2 ||V(A0 - F'm\ 

+c\\e\\\\Ae\\{\\v{A4> - F'[<t>))\\ + \\N>- F'm? + c) 

< ^WSuf + ^fMAt-F'mi 2 

+C(\\A<f> - F'(m 2 + IIA0 - F'(0)|| 4 + ||A0|| 2 + ||A#|| 4 ). 

The remaining terms can be estimated as follows 

K 3 <\a\\g\\\0\\\\Su\\<^\\Su\\ 2 + C\\A9f. 

lb 



K A < aA o ||V7r||||V0|| L 4||A</>-F'(0)|| L 4 
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< <7||Su||(||A0-F'(0)||5 +C)||A0-F , (0)||^||V(A0-F' 

< ^rll^H 2 + ^||V(A^ - F'(0))|| 2 + C(||A0 - F'(0)|| 2 + ||A</> - F"-" 1 - 
lb 8 



#5 ^aAoTll^'^llL-HA^-F'^)!! 2 < C||A</> - F'(0)|| 2 

K 6 = 2aA / Vi(A(j)- F'(<t>))VjUiS7j<l>dx 
Jn 

< C r ||V(A0-F , (^))[|||V U || I ,4||V0[| Il 4 

< ^||V(A0-F'(^))|| 2 + C||Vn||||An||(||A^-F'(0)|| + ||F'(0)||)||V0|| 

< ^ll^ll 2 + ^l|V(A<^ - F\m 2 + C(||Vn|| 2 + HV^II 4 + ||A0 - F'(</>)|| 4 ). 
lb 8 

It follows from the above estimates and the Sobolev embedding that 

^(HV^II 2 + aA o ||A0 - F'(^)|| 2 ) + f\\Suf + ^p||V(A* - F'(0))|j 2 

< C 5 ||VA0|| 2 + C(||Vu|| 4 + ||A0 - F'(</>)|| 4 + \\A6\\ 4 ) 

+C(||An|| 2 + \\A<j> - F(cJ))\\ 2 + ||A0|| 2 ). (3.4) 

We infer from (|1.15p that #t|r = 0, then it follows from (j!.13[) and the ^-equation ([l.lip that 
A8\r = 0. Applying A to both sides of (jl.lip . and taking the L 2 -inner product of the resultant 
with A9, we obtain 

- — \\A6\\ 2 + k\\VA6\\ 2 = -l A(u- V9)A9dx := K 7 , (3.5) 
2 ere Jn 



such that 



K 7 = f V(u-V6>)- VA6dx 
Jn 



< ^||VA0|| 2 + ^j|V(u- W)|| 2 

< ^||VA0|| 2 + C(||Vn|| 2 4||V0|| 2 4 + \\u\\loo || A0|| 2 ) 
8 

< -||VA0|[ 2 + C||Au||i||u||5(||A0||||W|| + ||W|| 2 ) + C||Au|||H|||Afl|| 2 
8 

< ^||VA0|| 2 + -^||5n|| 2 + C7|jVn|| 2 + C||A0|| 4 . (3.6) 
4 8KG5 



Hence, multiplying (|3.5p by 771 = kC§, and adding the result to (|3.4p . we deduce our conclusion 
()3.2p from ()3.6p . The proof is complete. □ 



Proof of Theorem SHI Since u G V, G H 2 (n),6 G (#o(0) n H 2 (Q)), we have 
.Ai(O) < +00. It follows from Corollary O that Ai{t) G 1^(0, +00). Then we infer from 
Lemma |3. II and [35^ Lemma 6.2.1] that Ai{t) is uniformly bounded for all time, which implies 

Kt)||Hi + ||^)l|^ + l|f(t)||tf»+ r +1 (ll^)ll^ + ll0(s)lli3 + ||0( S )||^)^<C, Vi>0. (3.7) 



where C is a positive constant depending on ||uo||hi> II^oIIh 2 ! II^oIIh 2 ) an d coefficients of the 
system. Then we can prove the existence of a global strong solution, which is actually unique 
by Proposition 13.21 The proof is complete. □ 
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For the weak solution, we still have Ai(t) € L 1 (0, +00). Then by ()3.2p and the uniform 
Gronwall lemma (cf. Temam [28} Lemma III. 1.1]), we conclude that for any 5 > 0, 

Ai(t + S) < c(i + ^\, Vt>0, 

where C is a positive constant depending on ||ito||, II^oIIh 1 ! II^oIIh 1 ! ^> an d coefficients of the 
system. As a result, 

Proposition 3.1 (Regularity of weak solutions in 2D). When n = 2, under the assumptions of 
Theorem 1 2. 11 any weak solution to problem (|1.8p - (|1.15p becomes a strong one for t > and the 
following estimate holds 

pt+i 

IK*)II^ + \\m\h + \W)\\b> + J (II«WIIh? + WMWh* + II'WIIhO* < v(t), vt > 0, 

where T> is a positive function depending on \\uq\\, H^oIIh 1 ) H^oH/p; an d coefficients of the 
system. In particular, lim T>{t) = +00. 

3.2 Three dimensional case 

Lemma 3.2. Suppose n = 3. Let (u,(j>,0) be a smooth solution to problem (|1.8p - (|1.15p . For 

the quantity A 2 {t) 

A 2 (t) = \\Vu(t)f + a\ \\Acf>(t) - F'i^tM 2 + \\&0(t)\\ 2 , (3-8) 
the following differential inequality holds 

j t A 2 (t) + u\\Su\\ 2 + aA O 7l|V(A0 - F'((f)))\\ 2 + A;||VA0|| 2 < C^A^t) + A 2 (t)), Vt > 0. (3.9) 

C* is a positive constant which may depend on \\uo\\, W^oWh 1 ) W^oWh 1 ! ^> an d coefficients of the 
system. 

Remark 3.1. We note that the coefficient of the third term in A 2 (t) is different from the one 
in Ai(t) (see (ETT]) ). 



Proof. We re-estimate the right-hand side of (|3.3p and (|3.5[) by using the 3D version of Sobolev 
embedding theorems. The estimates for K3 and K$ remain unchanged. Next, 

Ki < |M| l «>||Vu||||Sm|| < C\\Su\\i ||Vu|j§ < ^||5-u|[ 2 + C||Vn|| 6 . 



K 2a + K 2b < C\\Ve\\v\\Vct>\\l4Su\\ < — ||Su|| 2 + C||A0|| 2 ||A0|| 4 
< ^\\Su\\ 2 + C||A0|| 2 ||A^ - F'(0)|| 4 + C\\A9\\ 2 . 

K 2c < \b\X \\e\\ L o O \\A ( l>\\\\V(l>\\ Ij o a \\Su\\ 

< C , ||Ae||5||ve||3||A0||||V0||| ?l ||A0||5|| 1 S r «|| 

< Y^ll^ll 2 + C|| A0|||| V0||(||A0 - F\<P)\\ 3 + 1)(||V(A0 - F'(0))|| + ||A0 - F>(cf>)\\ + 1) 
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< ^ll^l| 2 + ^l|V(A0-F'(0))|| 2 

+C(||A0 - F'(0)|| 2 + \\A<f> - F'(</>)|| 8 + ||A0|| 2 + ||Afl|| 8 ). 

K A < a Ao||V7r||||V0|| L6 ||A^-F'(0)|| L 3 

< C\\Su\\(\\A0- F'(<P)\\ + 1) || - ^'(0)11 5 ||V(A0 - .P' 

< ^ll^|| 2 + ^p||V(A</> - F'(0))|| 2 + C(||A0 - F'(0)|| 2 + ||A0 - F' 

K 6 = 2aA [ V i (A0-F / (</.))n i V i V i (/)dx-2aAo(V(A</.-F / (0)),F / (0)n) 
Jn 

< C71|V(A0-i^(0))||(||ti|| L «||V 2 0|| L a + Hi^WllLooHull) 

< C||V(A</> - F'f»)||||Vu||(||A0 - F'(0)||i||V(A0 - F')\\* + \\A<f> - F>(<f>)\\ + l) 



< ^||V(A<£ - F'(0))|| 2 + C(||V U || 2 + ||Vnf + ||A0 - F>(4>)\\ 2 + ||A0 - F'(^)|| 8 ). 

O 

K 7 < ^||VA0|| 2 + i||V(u- W)|| 2 

< ^||VAfl|| 2 + C(||Vu|| 2 4 ||W|| 2 4 + H|l«||A0|| 2 ) 

< ^||VA0|| 2 + C||Su||i||Vu|[5||A0|| 2 



2 

*||VA0|| 2 + -, 
2 II II 12 



< -IIVA0P + — ||Sd| 2 + C||Vu|| 2 |A#ll 4 . 



Collecting all the estimates above, we arrive at our conclusion (j3.9|) . □ 

Proof of Theorem l3.2i Due to Lemma [3.2l a standard argument of the ordinary differential 
equation yields that there is a time To = Tq(u$, (f>o, 6o) < +oo such that Azit) is bounded 
on [0,Tq]. This enables us to prove that problem f)1.8p — fjl. 15j) admits a local strong solution. 
Uniqueness of the strong solution follows from Proposition 13.21 □ 

Since our problem contains the Navier-Stokes equation as a subsystem, in the 3D case, we 
cannot except the existence of global strong solutions to problem f)1.8p — fjl. 15j) for arbitrary large 
initial data. However, the global strong solution will exist if we further assume that the lower 
bound of the viscosity u* is sufficiently large. 

Set 

A 2 {t) = A 2 (t) + 1. 

Then Theorem 13.31 is a direct consequence of the following higher-order differential inequality 
concerning A 2 {t). 

Lemma 3.3. Suppose n = 3. We assume that v > 1 and (|2.2|) is fulfilled. Let (u,(p,9) be a 
smooth solution to the problem (|1.8p - (|1.15p . Then the following inequality holds 

j t M(t) +[u- M 1 ^A 2 (t)] \\Su\\ 2 + (aA o7 - Ml ^ 2 ® ) ||V(A0 - F'(0))|| 2 + k\\VA9\\ 2 
< M 2 A 2 (t), (3.10) 



Mi and M 2 are constants depending on \\uq\\, H^oIIh 1 ? H^oIIh 1 ! an< ^ coefficients of the system, 
but not on v . 



15 



Proof. We note that the uniform estimates in Corollary 12.11 still hold. Then we re-estimate the 
terms Ki,...,Ky in an alternative way. The estimates for K% and K§ remain unchanged. For 
the other terms, we have 

K x < C\\Su\\i\\u\\*\\Vu\\ 



< 



12 



\\Suf + z^||Vu|| 2 ||Au|| 2 + Cv-~\\Vu\\ 



^2a+^26<C||V0|| L 6||V0|| 2 6 ||^|| 

< C||VA0||^||V0||5(||A<A-F'(0)|| + l)(||V(A0-F'(0))||5||V0||Ul)!|5n|| 

< C||VAe||3||V6>||3([|A0-F'(0)|| + l)||5«|| 
+C||VA0||5||V0||5(||A0 - F'm + l)||Su||||V(A<£ - F'( 



< 



^ + uH\\M - F'(<p)\\ + ||W||)] \\Suf + ^||VA£|| 2 
+-(1 + ||A0 - F'(</>)|| 2 )||V(A<^ - F'(</>))|| 2 



+ c(l 



+ 



(||V(9|| 2 + ||A0-F' 



K 2c < |6|Ao||5||i=o||A0||||V0[|i,oo[|5«|| 

< C||A0||3||V0||3||A# (||VA0||3 + ||A0||l)||5n|| 

< C||A0||§||V0||§(||A0 - F'(<p)\\l + l)(\\V(A<p - F^))||3 + 1)||^|| 



< 



- + ^(||A0-F' 



+ ||A0||) ||5 U || 



+ (^ + i||A0-F^)|| 2 )||V(A0- J F' 
+C (\ + l)(||A0|| 2 + ||A</>-F' 



if 4 < C||^||(||A0-F'( ( A)|| + 1)||A0-F / (0)||2||V(A0-F'( 



< 



12 



+ ^||A0 - F'(<«|| 2 ) ||5 U || 2 + ^||V(A<£ - F'(0))|| 2 + -||A0 - F' 
/ 6 v 



Kq < C||A^-F'(0)|| L 3||Vu|| L 6||A^|| 

< C||V(A0 - F'((/)))\\^\\A<p - f (0)||*||Su||(||A0 - F' 

< ^ + ^\\A<p-Fi(cp)\\)\\Suf 

+ (^ + I||A0-F'(0)|| 2 ) ||V(A0-F' 



+C 1 + 



||A0-F' 



+ 1) 



K 7 



< 



< 



||VA#|| 2 + C(||V< 4 ||W|| 2 4 + 



MIL 



l|A0|| 2 ) 



||VA0|| 2 + C||^||2||Vn||2||A0|j2||V0||2 +C||5t*||5||«||2||Ae|| 5 
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< J||VAfl|| 2 + z^j|A0|| 2 ||Su|| 2 + -^-(||V-u|| 2 + ||A0|| 2 ). 



2 



Combining the above estimates, using the fact v > 1 and the definition of A 2 (t), we deduce the 
inequality (|3.10p . The proof is complete. □ 

Proof of Theorem 13.31 It follows from Corollary 12.11 that 

»t+l 

A 2 {r)dT < M, Vt > 0, (3.11) 

where M > may depend on Q, \\uq\\, \\(j)o\\ H i, I [ 1 1 j^ri ^ and coefficients of the system except v. 
Moreover, if the viscosity v is sufficiently large such that 

v\ > max jl, M i (A>(0) + M 2 M + 2M) + 1, (3.12) 

then following the same argument as in |17|.I31J , we can use Lemma 13.31 to obtain the uniform 
estimate 

~ , . z^2 min{aAn7, 1} 

A 2 (t)< ' W - ' (3 - 13) 

which yields the required conclusion. □ 

3.3 Uniqueness of strong solutions 

The uniqueness of strong solutions to the problem (|1.8|H|1.15|) can be obtained by the energy 
method. 

Proposition 3.2 (Uniqueness of strong solutions). For n = 2,3, let (m,0i,#i) and (u2,4>2,&2) 
be two strong solutions on [0, T] that start from the same initial data {uq, 4>o, 6q) £ V x H 2 (Q) x 
(-Hq(^) nH 2 (Q)) satisfying (fjE2j) . then (tii, </>i, 0i) = (u 2 , <fo, 2 ). 

Proof. We provide the proof for 3-D case and the proof for 2D is similar. Denote 

U = U\— 1i 2 , ^=01—02, 9 = 6\ — 6 2 - 

We can see that (u, ^, 0) satisfy 

(ut, u)v',V + / [C"i ' V)iti — («2 • V)u 2 ] • i>(ix + f / V-u : Vvdx 
Jn Jn 

= / [A(0i)V0i ® V0i - A(0 2 )V0 2 <g> V</> 2 ] : Vucte + ag / 0j • wcfcc, V« € V, (3.14) 
Jn ' Jn 

4>t + U!- Vfa - u 2 ■ V0 2 = 7 (A0 - F'(0i) + F'(0 2 )), a.e. in O, (3.15) 

t + ui • V0i - u 2 ■ V0 2 = kA6, a.e. in tt. (3.16) 

Taking v = u in d3HD , testing flHUgp by -aA o A0 and (|3~T6]l by -A(9 in L 2 (S1), respectively, 
adding up these three resultants, then performing integration by parts and using the incom- 
pressible condition for the velocity, we get 

--7-0NI 2 + aA o ||V0|| 2 + ||V0|| 2 ) + HlVull 2 + aAo7||A0|| 2 + £;||A0|| 2 
2 a I 

= -{u ■ Vui,u) - aA o (V<£A0i,u) + aA (ui • V4>,A<p) - 6A o (0V0i <g> V<pi,Vu) 
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-&A o (0 2 V0® V(f> u Vu) -6A o (02V02^V0,Vu) +ag(ej,u) 
+a\ o j(F'(0 1 ) - F'{<t> 2 ), A$) + (u • V0i, A0) + (u 2 • V0, A0) 
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:= Y, 1 ™- ( 3 - 17 ) 

m=l 

Keeping in mind the uniform estimates obtained in Corollary 12.11 we proceed to estimate the 
right hand side of (|3.17|) . 

h < NIl4||Vui|| < C||V«||f ||u||5||V«i|| < ^||Vn|| 2 + C||V-ui|| 4 ||u|| 2 . 

h < aA o ||V^|| L 6||iZ|| L 3||A0i|| < CIIA^IIIInll^lVnlllllA^II 
< ^||V^|| 2 + ^||A0|| 2 + C||A^|| 4 !|u|| 2 . 

h < |K|| L 6||A0||||V^|| L 3<C||V«i||||A0||I||V0||^ 
< ^||A^|| 2 + C||V Ul || 4 ||V0|| 2 . 



h + h < |6|Ao||Vu||||0|| L »||V</-i|| 2 4 + |a| 

< c||Vti||||ve|^||Ae||5 + c||v«||||V0|| 



< ^||VS|| 2 + ^||A^|| 2 + C||V0l 2 . 



h < |b|Aop2|U»||V«||||V0|| L 3||V0i|| L6 

< C||V«||||A^||5||V0||5||A0i|| 

< ^l|Vn|| 2 + ^||A^| 2 + C||A^|| 4 ||V0|| 2 . 



h < |6|Ao||fe|U»||V«||||V^||i I 3||V02||L« 

< a||Vu||||A^||i||V0||i||||A0 2 || 

< ^||Vn|| 2 + ^||A^|| 2 + C||A0 2 || 4 ||V0|| 2 . 



h < aA 7 



/ F"{8fa + (1 - s)<h)ds ||A0|| 



o 



< C||F"|MHM|A^|| 

< C||A0!|f||V0||^<^||A0|| 2 + C||V0|| 2 . 



h < ||«||i,3||Vei|| L 6||Ae||<C||Aei||||ti||i||V«||3||Ae|| 
< ^||V«|| 2 + ^||A0l 2 + C||A^|| 4 H| 2 . 



ho < ||« 2 || L 6||V0|| L 3||A^||<C||Vn 2 ||||V0||2||A^|| 

18 



< -||A£|| 2 + q|Vu 2 || 4 ||w|| 2 . 

6 



Summing up, we obtain that 



where 



j t {\\u\\ 2 + aAollV^I 2 + ||W|| 2 ) + v\\Vu\\ 2 + aA O 7l|A0|| 2 + k\\M\\ 2 
< Q(t)(||u|| 2 + aA o ||V0|| 2 + ||W|| 2 ), (3.18) 



Q(t) = C(l + ||V?/i|| 4 + IIA^H 4 + ||A^|| 4 + ||V« 2 || 4 + \\Acj) 2 \ 



with C being a constant that may depend on M, f2, and coefficients of the system. Since the 
V x H 2 x .ff 2 -norm of strong solutions to problem (|1.8|) - (|1.15|) are bounded on its existence time 
interval, Q(t) is bounded on [0, T]. Then we can complete the proof by applying the Gronwall 
inequality. □ 



Remark 3.2. (1) If n = 2, it is easy to check that f|3. 18|) holds with 

Q{t) = C(l + UVuill 2 + ||A0i|| 2 + ||A^|| 2 + ||Vu 2 || 2 + ||A</> 2 || 2 ). 

On the other hand, Corollarv \2.1\ implies that f+ + Q(s)ds < C for all t > 0. As a consequence, 
we can actually obtain the uniqueness of weak solutions to the problem (|1.8p - (|1.15p for the two 
dimensional case. 

(2) It is interesting to ask whether the problem (|1.8|) - (|1.15|) has a weak- strong uniqueness 
result in 3D as for the Navier-Stokes equations (c.f. Serrin \2J$ ). The main difficulty here is that 
due to the temperature- dependence of the surface tension coefficient X, we lose some dissipation 
in the derivation of the basic energy inequality fj2.3[) in order to control the corresponding higher- 
order nonlinear stress terms. When X is a constant, one can obtain the weak-strong uniqueness 
result as in fT 7 ^ for a simplified liquid crystal system. 



4 Long-time dynamics and stability 

In this section, we shall discuss the long-time behavior of global weak/strong solutions and 
stability properties of the problem (jl.8|) (|1.15|) . First, we present an alternative result that 
indicates the eventual regularity of weak solutions in 3D and is helpful to understand the long- 
time behavior of global solutions. 

Proposition 4.1. Suppose n = 3. For any (u o ,0o,#o) G V x H 2 (£l) x (Hq(£1) D H 2 (Q)) 
satisfying (|2.2p and 

MO) = llV^oll 2 + aA o ||A0 o - F'(0 O )|| 2 + j|A# || 2 < R, (4.1) 

where R > is a constant, there exists Eq > depending on \\uq\\, {{(PoWh 1 ) \\0qWh 1 ! ^> R> an d 
coefficients of the system such that either (i) the problem (|1.8|) - (|1.15|) admits a unique global 
strong solution that is uniformly bounded in time in V x H 2 x (Hq n H 2 ), or (ii) there is a 
G (0, +oo) such that £(%) < £(0) - e . 
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Proof. The proof follows from the idea in [17]. For the convenience of the readers, we sketch 
it here. Recalling the higher-order differential inequality (j3.9|) . we consider the following initial 
value problem of an ordinary differential equation: 

jY(t) = C*(Y 4 (t) + Y(t)), Y(0) = R> A 2 (0). (4.2) 
We denote by I = [0, T max ) the maximal existence interval of Y(t) such that lim Y(t) = +oo. 

ax 

It is easy to check that 

< A 2 (t) < Y(t), Vt € /, 

which indicates A(t) exists on /. We note that T max is determined by Y(0) = R and C* such 
that T max = T max (R, C* ) is increasing when R is decreasing. Taking to = \T max (R,C*) > 0, 
then we have 

< A 2 (t) < Y(t) <K, Vt € [0, t ], (4.3) 
where K is a constant that only depends on R, C*,io- Take 

e = ^Rt mm{v, 7, k}. (4.4) 

If (ii) is not true, namely, 6(t) > 6(0) — £0 for all t > 0, we infer from (|2.3p that 

y (2 ll Vu H 2 + aVy||A<£ - ^'(^ll 2 + k \\^ Q f) dt < e - 

Hence, there exists a t* G [§io>io] such that 

j / \ f2 1 li3e 

.A 2 (**) < max -,-,->—— = R. 

Iv 7 fc J to 

Taking £* as the initial time and Y(t*) = R in ()4.2p . then it follows from the above argument 
that Y(t) (and thus A 2 (t)) is uniformly bounded at least on [0, t* + to] ^ [0, §to]- Its bound 
remains the same as that on [0, to]- By iteration, it follows that A 2 (t) is uniformly bounded for 
t > 0. Thus, we can extend the (unique) local strong solution to infinity to get a global one. □ 

Proposition 4.2 (Eventual regularity of weak solutions in 3D). When n = 3, let (u,4>,8) be 
a global weak solution of the problem ()1 .8[ ) — f)l. 15[) . Then there exists a time Tq € (0, +00) such 
that (u,4>,9) becomes a strong solution in [Tq, +00). 

Proof. We simply take 

R = l, t = -^T max (R,C*), £ = -t min|^,7, fcj 



in the proof of Proposition 14.11 It follows from (|2.3p that there exist a T\ > such that 

"+00 

'Ti 



Then we can find a time T € [T x , Ti+§t ] such that A 2 (T ) < 1 and 6(t)-6(T ) > 6(t)-6(T{) > 
—£o for all t > To. Taking To as the initial time, we can apply Proposition 14.11 The proof is 
complete. □ 
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4.1 Convergence to equilibrium 

We shall show the convergence of global solutions to single steady states as time tends to 
infinity. Let (u, (j), 9) be a global weak solution of the problem (|1.8j) - (|1.15j) . We infer from either 
Proposition 13.11 (n = 2) or Proposition 14.21 (n = 3) that after a certain time T > 0, the weak 
solution will be a strong one that is uniformly bounded in V x H 2 x H 2 for all t >T. Since we 
are now considering the long-time behavior as t — > +oo, we can simply use a shift in time and 
reduce our study to the case of bounded strong solutions. 
The main result of this subsection is as follows: 

Theorem 4.1. Suppose n = 2,3. Any bounded global strong solution (u,(j),9) of the problem 
(|1.8p ~ (|1.15p converges to a steady state (0,</>oo,0) as time goes to infinity such that 

lim (\\u(t)\\m + \\cl ) (t)-^ oo \\ H 2 + \\9(t)\\ H 2) = 0, (4.5) 

t— >+oo 

where (poo satisfies the following nonlinear elliptic boundary value problem: 

- A0oo + F'Ooo) = 0, i£fi, urttfr 0oa| r = -l. (4.6) 
Moreover, we have the convergence rate 

K*)IIhi + U(t) -<ko||tf2 + \\Q(t)\\ H 2 < C(l + tyT^, Vt>0. (4.7) 

£ G (0,^) is a constant depending on 4>oo- Furthermore, 9 satisfies an exponential decay such 
that there exists a constant Co = Co(n, O) > ; 

< \\ d o\\e~ Cot , Vt>0. (4.8) 

Remark 4.1. Decay properties of the velocity u and temperature 9 as time tends to infinity can 
be obtained by the energy method (see Proposition ^ -3| below). However, convergence for the phase 
function <f> is usually nontrivial because the structure of the set of equilibria may be complicated 
and the solutions to elliptic problem like ()4.6p may form a continuum if the spatial dimension 
n > 2 (cf. e.g., Haraux fTH Remark 2.3.13]). As we have mentioned in the introduction, our 
results and their proofs hold for general Dirichlet boundary data for the phase function such 
that —1 in (|1. 14|) can be replaced by a certain generic function h{x). Since our problem enjoys 
a dissipative energy inequality (|2.3p . we can achieve the goal by using the Lojasiewicz-Simon 
approach (cf. e.g., fQ [T3\[T^[25\l ) . One advantage of this approach is that we can obtain the 
convergence result without investigating the structure of equilibria. 

The w-limit set of (tt , ^o,^o)GVx H 2 (VL) x (H^(Q) n H 2 (Q)) is defined as follows: 

u(uq, <p ,Oo) = {(u oo (x),0 oo (x),^ oo (x))GVx J ff 2 x( J ff 2 n J ff o 1 ): 
there exists {t n } oo such that 

(u(t n ),(j)(t n ),9(t n )) -»• (noo, 0oo,^oo ) in L 2 x H 1 x H 1 }. 

Proposition 4.3. Let n = 2,3. For any global strong solutions to problem (|1.8p - (|1.15p . there 
exists foo > such that 

lim £{t) = Soo, (4.9) 

t— »+oo 

and it holds 

lim (||«(t)|| H i + ||A0(t) - F'(m\\ + ||Afl(t)||) = 0. (4.10) 

t— >+oo 
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Proof. The total energy £(t) is nonnegative and decreasing as t increases (cf. (|2.3p ). Then 
(|4.9p easily follows. For global bounded strong solution (u,cj>,8), we have Ai(t) < C (n = 2) 
or A2(t) < C (n = 3), then it follow from Lemma 13.11 (n = 2) or Lemma 13.21 (n = 3) that 
< C (i = 1,2). On the other hand, we have known from ([273]) that G L^+oo). 

As a result, we can infer from Zheng |35[ Lemma 6.2.1] that limj_j. +00 Ai{t) = 0, which yields 
(14TTU1) . □ 

Corollary 4.1. cl>(moi </>0) $o) *s a nonempty bounded subset in V x iJ 2 x (.ff 2 fliJg). Moreover, 
u)(uo, 4>o, do) C S = {(0, 0,0)1 — Acf> + F'((f>) = in Q, (f>\r = — l} and f = ^ on u)(uq, 4>o, 6q). 

It is easy to verify that a critical point of the elastic energy E(<ft) given by (jl.7j) is equivalent 
to a solution to the following elliptic boundary value problem 



■ + F'((f))=0, ien, 0| r = — 1. (4.11) 
We recall the following Lojasiewicz-Simon type inequality (cf. Haraux-Jendoubi [13J). 

Lemma 4.1 (Lojasiewicz-Simon inequality). Let ^ be a critical point of E(cf>). Then there exist 
constants £ G (0, ^) and /3 > depending on ip such that for any cf) G i? 1 (Q) ; 0|r = —1 satisfying 

\\4> - iPWww < P> il holds 

|| - A<f> + F'(<l>)\\ H -i > \E(cj>) - E(t/i)\ l ~€, (4.12) 

For any global bounded strong solution, it follows from Corollary 14.11 that there is an in- 
creasing unbounded sequence {i n }ngN and a function (p^ G S such that 

lim Ufa) - tooWm = 0. (4.13) 
t„— H-oo 

As a result, we infer from Proposition 14.31 and (|4.13p that 

lim £{t n ) = aAo-E(^oo) = £oo and £{t) > aA -E(^oo), V t > 0. (4.14) 

t„— >+oo 

If £(i*) = aAo-E(</>oo) for some i* > 0, then £(t) = £oo for all t > 0. Thus, by (|2.3p . we have 
|[ n (t)|| v = ||0(t)|| H 2 = ||A0(t)-F'(^(i))|| = for t > t . Besides, it follows from equation (fTTTUD 
that 

Ut\\ < ||n|| L 4||V0|| L 4 + \\A<p - F'(cf>)\\ < C(\\Vu\\ + ||A0 - F'(<£)||), (4.15) 

thus ||</>t(i)|| = for t > t* and due to (|4.13p . we have cp(t) = 4>oo for t > to- 
Then we only have to consider the case £(t) > £oo = aX^E^^), for all t > 0. Based on 
the sequential convergence (|4.13p and the Lojasiewicz-Simon inequality, by using the classical 
argument in Jendoubi [14] . we can shown that after a certain time, the trajectory <f>(t) will fall 
into a certain small neighborhood of 0oo and stay there for all time. Namely, 

Proposition 4.4. There is a to > 0, such that \\<fi(t) — ^ooWh 1 < Pi f or a ^ t — to- 

Thus, for all t > to, <p(t) fulfills the condition in Lemma 14. 11 Since (1 - f) > \, then we infer 
from (gjgj) that 

(£(t) -aA o ^(0oo)) W < (\\u\\ 2 +(\\X79\\ 2 +u;\\e\\ 2 + a\o\E(<j))-E( 



■'oo) 



< C(\\u\\ 2 + ||W|| 2 + \\9\\ 2 ) l -t + C\E(<j>) - E^i 1 -^ 
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< C\\u\\ + C\\V6\\ + C||A0-F'( 
which combined with the energy inequality (|2.3|) yields that for t > to, it holds 

> ||Vu|| 2 + ||A0-F'(^)|| 2 + ||A^|| 2 



>|| + ||V0|| + ||A0-F'(0)|| H -i 
> C(||V«|| + ||A0-F'((^)|| + ||A0||). (4.16) 

Integrating (I4.16P with respect to t, using (|4.15p and the fact £{t) > £oo, we have 

/ Ut{t)\\dt < C / (||Vn|| + ||A<£ - F'(4>)\\)dt < C(£(t ) - < +oc, 

J to Jt 

which combined with the compactness of (j> in H 1 yields that lim^+oo \\4>(t) — (pooWii 1 = 0. 
Furthermore, since 

||A0- A0ool| < \\Acl)- A^-F' {(/>) + F'(^)\\ + \\F'((f>) -F'^W 

< ||A0-F'(</>)|| + ||<^)-<Ubi, (4.17) 



we conclude from (|4.10p that 



lim \\<p(t) - (pooWn 2 = 0. 

t— <r+oo 



It remains to prove the convergence rate (|4,7p . By Lemma 14.11 and f)4.16[) . we obtain that 

j t (£(t) - £00) + C{£(t) - £oo) 2(1 - C) < 0, Vt > t , 
which implies the decay rate for the total energy £ 

0<£{t)-£ oo <C(l + t)~ TJ ^, Vt> to- 
Integrating (I4.16P on (t, +00), where t > to, it follows from (|4.15p that 

r+OC e 

J Ut{r)\\dr < C(£(t) <C(l + t)-^3f. 
Adjusting the constant C properly, we get 

- 0oo|| < C(l + t)~^, Vt>0. (4.18) 



Higher-order estimates on the convergence rate can be achieved by constructing proper differ- 
ential inequalities via energy method. It is clear that for the asymptotic limit (0, <^oo,0), the 
system (|1.8p - (jl.lip is reduced to 

VPoo + ^VdV^ool 2 ) = -V0OO-A0OO, (4.19) 

-A0oo + F'^) = 0, with#oo|r = -l. (4.20) 

Denote tp = <j> — (p^. Then (u, if, 9) satisfies 

ut + u ■ Vit + Vp — vAu 
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= -aA o (Av9V0 + Ac^ooVv?) + 6A V • (6V4>® V</>) + a0gj, (4.21) 

V • u = 0, (4.22) 

+ « • = 7Ap - 7(^(0) - i^co)) , (4-23) 

(9 t + u • Vfl = £;A6>, (4.24) 

where we absorb all those gradient terms into the modified pressure p. 

Multiplying (jOTjl by u, g35D by aX (-Aip + (F'(<p) - F'(4> 00 ))) + if and (Oil by -A0 
respectively, integrating over Q and adding them together, we have 

\jt m + H|Vn|12 + aX °^ " F '^H 2 + + k W 

= -b\ 9V i (j>Vj<l)VjUidx + ag(6j,u) + (u-'V6,A6) 
Jn 

5 

:= ^i? m , (4.25) 

m=l 

where 

3>(i) = INI 2 + aAo||V<^|| 2 + 2aA [ [F(<f>) - F^) - F ' (M^dx + |M| 2 + ||V0|| 2 . (4.26) 

Jn 

In the derivation of (|4.25p . we have used (|4.19|) . (|4.2U|) and the following fact 

/ (AipVcp + A^ooVc^) • udx + / u-V(j)[-Aip+ (F' ((/))- F'(^oo))]dx 
Jn Jn 



/ (Ac^oo - F'(<ko)) V0 • udx + / u ■ VF(^)dx - / A^V^ 
Jn Jn Jn 



udx 



0. 



Since we are now dealing with global strong solutions that are uniformly bounded in V x H x H , 
it follows that 

Ri < ||^||^||V U ||||V0|| 2 4 <C||A0||l||^||i||Vn||||^||^ 

< ^\\Vuf + t\\A6\\ 2 + C\\ef, 

R 2 < \a\\g\\\6\\\\u\\<^\\Vu\\ 2 + C\\9\\\ 

R, < ^||A^|| 2 + C|| U || 2 6 ||V^|| 2 3<^||A0|| 2 + C7P||5||A0||I 

< ^||A0|| 2 + C||0|| 2 , 

«4 + «5 < Ikll^llV^llLall^ll+Cll^f^^llVnf + CI^H 2 . 

From the definition of F(<p), we have — -^(^oo) — F'((j) 00 )ip]dx\ < Cill^ll 2 - Combined 

with the definition of y(t) in (|4.26|) . it yields that 

y(t) + d|M| 2 > C(\\u\\ 2 + + ||W|| 2 ). (4.27) 

It follows from ()4.25p . ()4.27p and the estimates on R m (m = 1, 5) that 

j t y(t) + c 2 y(t) + c 3 Mt) < c(y\\ 2 + \\ef), % = 1, 2. (4.28) 



24 



Recalling Lemma 13. II (n = 2) or Lemma 13.21 (n = 3), we have 

d 



di 



< C 4 Ai{t). 



Multiplying (|4.29|) with r\ = l^h , and adding the resultant to f|4.28j) , we get from 



(4.29) 
and flUBD 



that 



dt 



[y(t) + T)Mt)] + c'[y(t) + T)Mt)] < c(IMI 2 + ll^ll 2 ). v * > °- 



Consequently, 



y(t) + v Mt) 



< Ce- C,t + Ce- C,t 



Ce-c't + Ce- '* 



e C ' T (l + r) i-29 dr + / e°" r (l + r) i-aerfr 



e c ' T (l + T)-^sdT 



+Ce 



-Ct 



,G't 



c 



-(1 + t)' 



29 
' 1-29 



T=t 



+ 



29 



C'(l - 29) J 



e c " T (l + r) i-29 dr 



2€ 



< c{i + t)-—t, 



( - I (1 + r) 1-29 dr + (1 + i) i-29 e c 
o 

Vt > 0. 



(4.30) 



Then our conclusion (|4.7p follows from (|4.30p . the definitions of Ai(t), y(t) and (|4.17p . The 
exponential decay of 9 (see (|4.8p ) easily follows from (|2.6p . The proof of Theorem l4.1l is complete. 



4.2 Stability of energy minimizers 



We have shown that any global weak (or strong) solution of problem (|1.8p - (|1.15|) will converge 
to a steady state as time goes to infinity (without smallness restrictions on the initial data). 
However, it is not clear to which equilibrium it will converge, since the set of equilibria may be 
a continuum (for general Dirichlet boundary data of the phase function) . This is different from 
the classical concept of stability in the literature. Below we shall show that if no and 9q are close 
to zero and </>o is near a certain local minimizer of the elastic energy E(<p), problem ([1.8p - ([1.15p 
admits a unique global strong solution. Moreover, the energy minimizer is Lyapunov stable. 
The same as in Remark 14.11 the results and their proofs in this subsection actually hold for 
general Dirichlet boundary conditions for (p, not only the special case (|1.14p . 

Definition 4.1. The function eft* G /C := {(f) € H l (Q) : <^>|r = —1} is called a local minimizer 
of E((p), if there exists a > 0, such that for any (j) € /C satisfying ||0 — </>*||#i < o, it holds 
E(4>) > E{4>*). 

Remark 4.2. It is easy to verify that any local minimizer of E((p) is a critical point of E((p) 
and satisfies the elliptic boundary value problem (|4.1ip . 

Theorem 4.2. Suppose n = 3 and (|2.2p is satisfied. Let 4>* G H 2 (£l) n K, be a local minimizer 
of E{4>). For arbitrary r > 0, we consider the set 



B r = {(u,4>,6) G V x (H 2 (n) n/C) x (H 2 (n) 0^(0))} : 
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|M| H i < r, U ~ <t>*\\H* < r, \\e\\ H 2 < r}. (4.31) 

Then there exist positive constants ci,^,^ depending onr, Q, a, v* , <f>* and coefficients of the 
system, such that for any initial data {uq, 0o> &o) € &r satisfying 

\\uq\\ < ai, Ho - 4>*\\m < <?2, Po\\m < 0-3, (4.32) 

we have 

(i) the problem (|1.8p - (|1.15p admits a unique global strong solution (u,(f>,9); 

(ii) the energy minimizer <p* is Lyapunov stable; 

(Hi) the global strong solution has the same long-time behavior as in Theorem \4.1\ Although 
the limit function (p^ may differ from the minimizer (j)* , the total energy £{t) will converge to 
the same energy level of (f>* 

lim £(t) = aA o £(0oo) = aA £(</>*). (4.33) 

t— »+oo 

Moreover, if <fi* is an isolated local minimizer, then (f)^ = (j)* and thus (j)* is asymptotically 
stable. 

Proof. By Proposition 14.11 in order to prove the existence of global strong solutions, we only 
have to verify that 

S(t) -£{0) > -eo, Vt€[0,+oo), (4.34) 

where eo is defined as in (|4.4p . We know from the argument in Proposition 14.11 that there exists 
to = jT max and A2(t) is uniformly bounded on [0, to] by a constant depending on r, <f>*, U and 
coefficients of the system. Since Azft) is bounded on [0,to], it holds 

£(t)-£(0) 

= \\Wt)\\ 2 - \\\u \\ 2 + a\ (Em)) ~ E{<M) + ^||W(t)|| 2 + |||0(t)|| 2 

-§l|V<y 2 -f|W| 2 

> ~\\u \\ 2 + aA - E(<f>*) + E{<jf) - E{fa)) - |||V0 O || 2 - ||N| 2 

> -^INII 2 - |||V0 O || 2 - ||N| 2 ~ C x \\<h - <t>*\\ H i + aA o (£(0(t)) - E(f)), 

where Q,uj are as in Proposition 12. II and C\ depends on r, <f>* , f2, and coefficients of the system. 
We assume that a m (m = 1, 2, 3) are sufficiently small such that 

-a 2 + - max{C, w}cr| + Cia 2 < £q. (4.35) 

If we can ensure that 

E{<t>{t)) - E(<f>*) > 0, V t e [0, t Q ], (4.36) 

then we have 

£(t)-£{0) > -e , VtG [0,t ]. (4.37) 

This enables us to apply the argument in the proof for Proposition 14. 1 1 to extend the local strong 

fL] = [0, f> 



solution from [0, t ] to [0, t + ^] = [0, ^1 
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By Definition 14.11 (|4.36p can be reduced to the following condition 

\\<f>(t)-<l>*\\ m <mm{<T,p}:=5, VtE[0,t ], (4.38) 

where f3 > is the constant depending on (jf given in Lemma [4.11 (taking ip = <ft* therein, we note 
that <f>* is a critical point of E(<p)). We shall show that one can choose a smaller 02 satisfying 

0-2 < -6, (4.39) 

such that ()4.38p holds. This can be done via the Lojasiewicz-Simon approach by a contradiction 
argument (cf. Wu et al [32J). If (|4.38p is not true, then by the continuity of (j> that <f> € 
C([0, to]; Hq), there exists a minimal time To G (0, to], such that ||0(To) — 0*||#i = 5. We 
observe that £(t) > a XoE ((/>*) for any t € [0, To]. If for some T < To, £(T) = a\oE((j)*), then we 
deduce from the definition of the local minimizer and the basic energy inequality (|2.3p that for 
t > T, £(t) cannot drop and will remain a\ E((f)*). Thus, Vu = A<j) - F'(<f>) = A6 = for all 
t > T and the evolution becomes stationary. The conclusion easily follows. In the following, we 
just assume £(t) > a\$E((f)*) for t € [0, To]. Applying Lemma UTTI with ip = <f>* , we get (similar 

to (jUSD ) 

- ^[5(t) - a\ E{<t>*)f > C(\\Vu\\ + ||A0 - F'(<f>)\\ + \\A0\\), Vt G (0,T ). 
Then we infer from (|4,15p that 

U(T ) - Mm < C||<KTq)-0o||%(T o )-0o||| 2 

< c(jf <C[f(0)-oAoW)]§ 

< C 2 (\\u \\ 2 + \\6 \\ 2 H1 + ||0 O -^|| H1 ) 5 . (4.40) 
Choosing <r m (m = 1,2,3) satisfying ([Q5]> . (|4T39l) and 

C 2 (^ + a 3 2 + a 2 ) l <^, (4.41) 

we can see that 

||0(T O ) - 0*11^1 < ||0(T O ) - ^qIIhi + ll^o - <\s<S, 

which leads to a contradiction with the definition of Tq. Thus, ()4.38p is true and (|4.37p holds. 

By iteration, we conclude that the local strong solution (it, cf>, 9) can be extended by a fixed 
length ^ in each step and it is indeed a global solution with A^ify being uniformly bounded. 
Then by Theorem I4.lt there exists a critical point (j)^ of E(cp), such that the global solution 
(u(t), 4>(t),6(t)) satisfies the same long-time behavior (|4.5p with convergence rate (|4.7p . It is 
easy to see from the above argument that for any e > 0, by choosing sufficiently small a m 
(m = 1,2,3), it holds \\4>{t) — 4>*\\hi < e, for all t > 0. This implies the Lyapunov stability of 
the local minimizer <j>*. In particular, we have 

Woo - (f>*\\m < \\4>{t) - ^Whi + \\<t>(t) - 4>*\\ m < min{<7,/3}. 

Applying the Lojasiewicz-Simon inequality once more with ip = </>*, we conclude that 

1^(000) - £(0*)| W < || - + F'i^W = 0, (4.42) 

which together with (|4.5p yields ()4.33p . The proof is complete. □ 
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